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ABSTRACT 

To measure the column densities of interstellar and interg alactic gas clouds using a bsorption line 
spectroscopy, the apparent optical depth technique (AOD) of iSavage fc Sembachl (|10Qlf ) can be used 
instead of a curve-of-growth analysis or profile fit. We show that the AOD technique, whilst an excel- 
lent tool when applied to data with good S/N, will likely overestimate the true column densities when 
applied to data with low S/N. This overestimation results from the non- linear relationship between 
the flux falling on a given detector pixel and the apparent optical depth in that pixel. We use Monte 
Carlo techniques to investigate the amplitude of this overestimation when working with data from 
the Far Ultraviolet Spectroscopic Explorer (FUSE) and the Space Telescope Imaging Spectrograph 
(STIS), for a range of values of S/N, line depth, line width, and rebinning. AOD measurements of 
optimally sampled, resolved lines are accurate to within 10% for FUSE /LiF and STIS/E140M data 
with S/N>7 per resolution element. 

Subject headings: techniques: spectroscopic 
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1. INTRODUCTION 

The apparent o ptical depth (AOD) method 
ijSavaee fc Sembacbl Il991t iSembach fc Savaed 11992(1 
has gained substantial popularity as a means of convert- 
ing velocity-resolved flux profiles into column density 
measurements for interstellar and intergalactic absorp- 
tion lines. The method offers a quick and convenient 
way of determining reliable interstellar column densities 
from unsaturated lines without having to follow a full 
curve-of-growth analysis or detailed component fit, and 
without demanding prior knowledge of the component 
structure. Accurate column densities are important for 
studies of elemental abundances and physical conditions 
in diffuse gas in space. 

In the AOD method, a velocity-resolved flux profile 
F(v) is converted to an apparent optical depth profile 
T a (v) using the relation: 

Ta (v)=\n[F c (v)/F(v)], (1) 

where and F(v) and F c (v) are the observed line and con- 
tinuum fluxes at velocity v, respectively. This apparent 
optical depth profile can then be converted to an appar- 
ent column density profile according to 

N a (v) = 3^768 Xl0 14 (/A) -1 T o (i>) cm -2 (kms^ 1 )- 1 , (2) 

where / is the oscillator strength of the transition 
and A is the transitio n wavelength in Angstroms 
l)Savage fc Sembacbl 1 19911) . The total apparent column 
density between two velocity limits v— and v+ is then 
simply N a = J v _ N a (v)dv. When two lines of different 
strength of the same ionic species are available, the AOD 
method can also be used to assess and correct for the 
level of saturation in the data, by comparing the appar- 
ent column density derived from th e stronger line with 
that derived from the weaker line (Savage fc Sembacbl 

The AOD method was originally developed for applica- 
tion to measurements of reasonably high signal-to-noise 
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spectra (S/N>20 per resolution element). The technique 
offers excellent results in these cases. However, spectro- 
scopists are now using the m ethod even when an alyzing 
relatively low S/N data (e.g. lWakker et alJl20^ . This 
paper deals with accounting for a systematic error that 
is introduced when applying the AOD method to noisy 
data, which can lead to an overestimation of the true col- 
umn density of the absorbing species. The error arises 
because the logarithmic relationship between T a {v) and 
F(v) (Eq. 1) will distort a Poisson noise component in 
the flux when converting into column density space, giv- 
ing undue weight to the pixels where the noise has re- 
sulted in high optical depths; this effect tends to exag- 
gerate the estimate of N a . We present simulations inves- 
tigating the accuracy of N a determinations as a function 
of line depth and S/N, together with the separate effects 
of line width and rebinning. 

2. SIMULATIONS 

We ran Monte Carlo simulations of the measurement of 
noisy absorption lines to investigate the amount by which 
noise in the data leads to overestimation of N a . Our 
models were run with two simulated experimental setups, 
representing t he Far Ultraviolet Spectroscopic Explorer 
(FUSE) LiF ijSahnow et alJl200C|) and Space Telescope 
Imag ing Spectrograph (STIS) E140M fWoodgate et alJ 
1998) configurations. The FUSE /LiF setup is parame- 
terized by a Gaussian instrumental line spread function 
with FWHM=20kms~ 1 (a ins = S.^kms^ 1 ) and veloc- 
ity pixels 2.0kms~ 1 wide. The STIS/E140M configu- 
ration has FWHM^e.Skms" 1 (<r lns = 2.9kms" 1 ) and 
3.2 km s -1 pixels. Our models simulate the measurement 
of an absorption line whose intrinsic optical depth profile 
is a single-component Gaussian, centered for convenience 
at Okms" 1 , and normalized to a peak optical depth of 



to, i-e. 



t(v) 



r e' 



~ v / 2<T ii,, e . We run models for two 
ine = 2cri ns ) and marginally 



line width cases: resolved (ai 

resolved (au ne = Cjns) The AOD method is known to 
underestimate the column density for unresolved lines, 
when no allowance is made for the effects of unresolved 
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saturation [see Table 5 of iSavage fc Sembachl 1(1991')]. 

For a grid of different intrinsic central line depths du ne 
we converted the optical depth profile to a normalized 
flux profile using F norrn = e~ T ^ v \ where the intrinsic 
central optical depth To = — ln(l — du ne ). The flux pro- 
files were then convolved with the instrumental broad- 
ening function; the effect of this process is to decrease 
the intrinsic line depth to an apparent line depth (<i), 
and increase the intrinsic line width to an apparent line 
width (a). We then added different levels of random 
noise (drawn from a Gaussian distibution) to the data, 
scaled so that the S/N in the continuum reached the 
desired levels when rebinned to pixels equal in size to 
the resolution element. In the line, we scaled down the 
random noise contribution by \/F, in accordance with 
Poisson statistics. 

We then rebinned the data by different factors, con- 
verted the profiles to N a (v) space, and finally integrated 
the profiles over velocity to yield N a , as a function of 
line depth, noise, and rebinning factor. The velocity in- 
tegration limits used are ±2.5<r around the line center, 
corresponding to the points where the line recovers to 
96% of the continuum in the absence of noise. These 
limits were chosen so as to reproduce the velocity inte- 
gration ranges used in practice, determined by compar- 
ing the line widths and velocity integration ranges from 
a large spectroscop ic data set (the FUSE O VI survey; 
Wak ker et alJEOOfl) . For each model run, we determined 
the ratio of N a to N trU e, where N trU e is the true col- 
umn density of the simulated component, then repeated 
the process 500 times; the mean overestimation factor 
< N a > /Ntme, together with its dispersion, is then 
used in our results. 

For very optically thick lines, the random error can 
take the flux below zero, corresponding to negative 
N a (v), which has no physical meaning. Therefore, a cut- 
off must be used, whereby all points with F(v) less than 
some value are placed at that value (in our case, 1% of 
the continuum). The cutoff value selected will depend 
on the reliability of the scattered light and background 
corrections for the particular instrument providing the 
observations. 

3. RESULTS 

Figure 1 shows a graphical example of the overestima- 
tion process for a resolved line with an intrinsic depth of 
0.8, measured with FUSE. In the left column we show 
how the spectra would appear at different S/N levels, 
ranging from S/N=oo at the bottom to S/N=5 at the 
top. Note that we quote the S/N per resolution element, 
rather than the S/N at the level of rebinning shown. The 
two are related by (S/N) res e ; em .=^/(number of bins per 
resolution element) x(S/N)(,j„. We verified this relation 
by measuring the S/N of the model spectra at different 
binning levels. In the next three columns we show the 
profiles in apparent column density space after different 
levels of rebinning. 

In Figure 2 we display the results of different pa- 
rameter cases using a color-coded plot: resolved lines 
with FUSE (top left), marginally resolved lines with 
FUSE (top right), resolved lines with STIS (bottom left), 
and marginally resolved lines with STIS (bottom right). 
Within each panel, we show the overestimation factor as 
a function of signal-to-noise ratio and apparent line depth 



(both easily measurable quantities), for three different 
rebinning cases. In each panel, we have not shaded in re- 
gions that correspond to non-significant detections, i.e. 
those with W\ < 3a(W\); the minimum depth line that 
can be considered significant is proportional to (S/N) -1 , 
since at low S/N shallow lines will not be distinguished 
from the noise. Using Figure 2 an observer can find the 
likely error on N a in a line with a given depth, width, and 
S/N, and correct for it. The overestimation factors for a 
range of cases are also presented in Tables 1 and 2, for 
FUSE/UF and STIS/E140M data, respectively. Upon 
request we can also provide an IDL code that returns 
the overestimation factor when given the measured S/N, 
apparent line depth, and apparent line width. 

4. DISCUSSION 

The effects of several parameters combine to produce 
the trends seen in Tables 1 and 2 and Figures 1 and 2: 
S/N, level of rebinning, spectral resolution, line width, 
and line depth. We discuss each separately in this sec- 
tion. 

• Low S/N and Rebinning - the non-linear dis- 
tortion of the N a (v) profile and resulting overesti- 
mate of N a becomes worse as the S/N decreases 
(Fig. 1). Rebinning the spectrum before measure- 
ment can lessen the distortion, since the rebinning 
process increases the effective S/N. Once the data 
is rebinned to a level coarser than the resolution el- 
ement, the gain in the N a accuracy is made at the 
expense of resolution, so that kinematic informa- 
tion is given up. However, if the absorption line be- 
fore rebinning is fully or marginally resolved, and if 
one is more interested in obtaining an accurate col- 
umn density than obtaining knowledge about the 
detailed shape of the line profile, then overbinning 
is justified at low S/N. AOD measurements made 
at optimum sampling (two rebinned pixels per res- 
olution element) on FUSE and STIS data are over- 
estimated by less than 10% for S/N per resolution 
element >7. For data with lower S/N, rebinning 
beyond the resolution element will likely increase 
the accuracy of the derived column density. 

• Lack of Resolution and Line Width - for 
marginally resolved lines, the apparent column 
density tends to underestimate the true column 
density at high line depths (purple shading on the 
right hand panels in Figure 2), since the high op- 
tical depth points are smoothed out to lower ap- 
parent optical depths by the instrumental blur- 
ring. This effect therefore works in the opposite 
sense to the S/N effect. The more unresolved 
the line is (i.e., the lower crime/ a ins), the worse 
this effect becomes. Even full saturation can be 
hidden by the effects of i nstru mental broade ning. 
ISavage fc Sembachl l)1991|) and I Jenkins! (|1996|) dis- 
cuss how this unresolved saturation can be cor- 
rected for when observations exist for more than 
one line of an ion. For resolved lines, increasing 
the ratio cru ne /ai na above 2 increases the number 
of pixels susceptible to noise and hence increases 
< N a > I N true slightly for a given S/N and line 
depth (by a few percent if crime/ Oins — 4); how- 
ever, since the dependency on width is weak, we 
present results for two line width cases only. 
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• Line Depth - for a given S /N and line width, the 
overestimation factor is worse for shallow lines than 
for deep lines. This is because the noise is Poisso- 
nian, so a pixel with low optical depth (shallow) 
has a higher error on the flux than a pixel with 
high optical depth (deep). For very deep lines, the 
cutoff used to eliminate negative pixels underesti- 
mates the optical depth in very opaque parts of the 
profile. This effect works in the opposite sense to 
the S/N effect, so that measurements of very deep 
(F(v)/F c (v) < O.fO) lines at low S/N are distorted 
by two competing effects. Note that if the noise 
were random rather than Poissonian (i.e. indepen- 
dent of the flux), the overestimation effect on deep 
lines would become worse since the flux errors deep 
in the lines are larger. 

In Figure 3 we show the histogram of N a /N trU e for 
one particular run (FUSE/hiF, <Ju ne — 2cri ns , S/N=8, 
intrinsic depth=0.5, 500 runs) for three levels of rebin- 
ning. The overestimation factor has a Gaussian distri- 
bution, with a mean that decreases with rebinning. For 
any given absorption profile, the overestimation can take 
a range of values. We therefore advise a conservative 
error estimate when correcting for the overestimation. 



We note that the non- linear S/N distortion effect dis- 
cussed here can influence the detailed intercomparisons 
of N a (v) profiles of the same ion. Therefore, corrections 
for line saturation based on comparing N a (v) profiles of 
weak and strong lines may not be valid for low S/N ob- 
servations. 

In their survey of O VI absorption in the Ga lactic halo 
and high-velocity clouds, iWakker et alJ l)2003j) rebinned 
FUSE data with S/N<10 per resolution element by 5 
pixels, and noisier data by up to 20 pixels, before mea- 
surement with the AOD technique. Our modelling shows 
that this was the correct procedure for measuring accu- 
rate O VI column densities, and we endorse a similar 
treatment of low S/N data in future studies. Without 
this rebinning, the column density measurements along 
sight lines with noisy spectra could be inflated by more 
than 20%. 
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TABLE 1 
Simulation Results - FUSE /LiF 



Intrinsic 




Apparent 




S/N b 




< Na/Ntrue > 




FWHM 


Depth 


FWHM 


Depth 




3 pixel rebin 


5 pixel rebin 


10 pixel rebin 


(kms" 1 ) 


(%) 


(kms- 1 ) 


(%) 




(6 km s~ 1 bins) 


(10kms~ 1 bins) 


(20 kms" 1 bins) 


/in n 
4U.U 


on 


A K Q 
40. o 


1 Q 
lo 


i n 
1U 


l.zuztu.01 


111 4~n A o 


1 0.74-0. tro 










8 


1.32±0.o9 


l.lo±0.64 


1.09±0.o8 










6 


1.73±1.05 


1.37±0.91 


1 O A I O 1 

1.24±(J.91 










4 


O 1 O I O 1 ,4 


o oo I i 7n 


1 A 0_L_ 1 A 7 

1.49±1.47 


,in n 
40.0 


50 


4/. 5 


45 


i n 
10 


1 fi7±n on 
l.U / ±U.2U 


1 n /i _i_n i n 

l.U4±U.ly 


1 01 _lo on 
l.Ul±U.2U 










8 


1 in 1 n o T 

1.13±0.27 


1 OCT I O O A 

1.05±0.24 


1 1 I O O A 

l.Ul±U.z4 










6 


1.29±0.44 


1 1 1 I n OO 

1.13±0.38 


1 f\ A I n or 

1.04±0.35 










4 


1.92±0.77 


1 A n I o 7n 

1.49±0.70 


1 on I O ET£? 

1.20±0.5o 


40.0 


99 


61.1 


97 


10 


0.87±0.U7 


0.85±0.U7 


n on I n no 

0.80±0.U8 










8 


o on I n nn 

0.90±0.U9 


n o T I n nn 

0.87±0.09 


n o 1 I n 1 n 

0.81±0.10 










6 


o oo I o 1 n 

0.93±0.10 


n on I o 1 i 

0.89±0.11 


n o o I n i o 

0.83±0.12 










4 


i nn I n 1 r* 

l.UOiO.lb 


0.94±0.1o 


n o rr I n 1 T 

0.85±0.17 


20.0 


20 


28.6 


n 

14 


10 


"I "I I t'\ \. > "I 

1.21±0.81 


1.11±0.81 


1.04±0.87 










8 


1.35±1.03 


1.19±1.01 


1.05±1.06 










6 


1.78±1.57 


1.41±1.44 


1.18±1.44 










4 


3.42±3.35 


2.15±2.71 


1.45±2.34 


20.0 


50 


29.3 


37 


10 


1.03±0.30 


0.99±0.31 


0.94±0.30 










8 


1.11±0.42 


1.04±0.41 


0.97±0.40 










6 


1.27±0.62 


1.12±0.56 


0.99±0.52 










4 


2.06±1.34 


1.56±1.18 


1.18±0.92 


20.0 


99 


34.2 


89 


10 


0.67±0.09 


0.63±0.08 


0.56±0.07 










8 


0.69±0.12 


0.65±0.11 


0.57±0.08 










6 


0.75±0.16 


0.69±0.15 


0.58±0.12 










4 


0.88±0.25 


0.79±0.26 


0.64±0.23 



Note. — The FUSE /LiF spectrograph has FWHM;„ 3 =20 km s 1 and pixel size 2.0 kms 1 . The upper and lower halves of this table represent 
fullly resolved and marginally resolved lines, respectively. 

instrumental broadening causes the apparent line depth to be less than the intrinsic line depth, and the apparent width to be broader than the 
intrinsic line width. 
b S/N is quoted per 20 kms" 1 resolution element. 
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TABLE 2 
Simulation Results - STIS/E140M 



Intrinsic 




Apparent 


a 


S/N b 




< Na/Ntrue > 




FWHM 


Depth 


FWHM 


Depth 




1 pixel rcbin 


2 pixel rcbin 


3 pixel rcbin 


^km s ^~ ) 




(kms -1 ) 








( (\ A. Lrtn t; ^ ViintA 




13.6 


20 


15.6 


18 


10 


1.11±0.52 


1.04±0.54 


1.03±0.53 










8 


1.16±0.70 


1.04±0.70 


1.02±0.71 










6 


1.48±1.00 


1.28±0.97 


1.21±0.94 










4 


2.11±1.81 


1.46±1.52 


1.29±1.43 


13.6 


50 


16.2 


45 


10 


1.03±0.18 


1.00±0.19 


0.98±0.19 










8 


1.06±0.25 


1.01±0.24 


0.98±0.23 










6 


1.11±0.36 


1.02±0.34 


0.99±0.33 










4 


1.44±0.69 


1.14±0.53 


1.06±0.48 


13.6 


99 


20.8 


97 


10 


0.87±0.08 


0.82±0.08 


0.75±0.07 










8 


0.89±0.09 


0.84±0.09 


0.76±0.09 










6 


0.91±0.11 


0.85±0.12 


0.78±0.13 










4 


0.96±0.16 


0.87±0.17 


0.81±0.18 


6.8 


20 


9.7 


14 


10 


1.10±0.78 


1.02±0.81 


1.01±0.85 










8 


1.21±1.05 


1.09±1.08 


1.05±1.16 










6 


1.43±1.45 


1.24±1.48 


1.17±1.56 










4 


2.50±2.98 


1.76±2.48 


1.55±2.35 


6.8 


50 


10.0 


37 


10 


0.98±0.28 


0.94±0.30 


0.92±0.31 










8 


1.04±0.37 


0.99±0.37 


0.95±0.38 










6 


1.09±0.51 


0.98±0.49 


0.93±0.51 










4 


1.52±1.10 


1.18±0.87 


1.09±0.86 


6.8 


99 


11.6 


90 


10 


0.66±0.09 


0.60±0.07 


0.53±0.07 










8 


0.67±0.12 


0.60±0.10 


0.52±0.09 










6 


0.72±0.16 


0.62±0.14 


0.54±0.12 










4 


0.84±0.26 


0.69±0.25 


0.58±0.21 



Note. — STIS/E140M has FWHMi„,=6.8kms _1 and pixel size 3.2 km s . The upper and lower halves of this table represent fullly resolved 
and marginally resolved lines, respectively. 

instrumental broadening causes the apparent line depth to be less than the intrinsic line depth, and the apparent width to be broader than the 
intrinsic line width. 
b S/N is quoted per 6.8 kms -1 resolution element. 
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Spectra: rebin by 3 



- S/N 



a 



Na(v): rebin by 3 
Oversampled 





Na(v): rebin by 5 
Optimum Sampling 



-100 100 
Velocity 



— 0- 
q.94±0 i .00 

-50 50 
Velocity 




Na(v): rebin bylO 
Undersampled 



— & 
q.93±0,00 

-50 50 
Velocity 




0.92±0.00 

-50 50 
Velocity 



Fig. 1. — Simulated spectra illustrating the effect by which N a is overestimated when the AOD technique is applied to noisy data, for 
resolved absorption lines observed with FUSE with an intrinsic central depth of 80%. The left column shows the simulated spectra with 
various levels of Poisson noise added. The next three columns show the effect of translating the noise into apparent column density space, 
after different levels of rebinning. The N a (v) profiles have been arbitrarily normalized for convenience. The numbers written in the corner 
of each panel represent the mean column density overestimation factor < 7V a > /Ntrue after 500 runs, together with its dispersion, where 
N a = / N a (v)dv and Ntrue is the true column density in the absorber. S/N is quoted per resolution element (corresponding in this case 
to the S/N at 10 pixel rebinning). 
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Fig. 2. — Color-coded plots showing the column density overestimation factor N a /Ntrue as a function of signal-to-noise per resolution 
element, apparent line depth, and rebinning, for absorption lines measured using the AOD method in four cases: FUSE resolved lines 
(top-left), FUSE marginally resolved lines (top-right), STIS resolved lines (bottom- left), and STIS marginally resolved lines (bottom-right). 
Only lines with > 3<r equivalent width detections are included. Green regions show cases where the AOD method returns the column density 
correct to within 20%; blue and purple regions show where the AOD method will underestimate the true column density (unresolved lines, 
particularly at high depth); yellow and red regions show cases where AOD measurements lines will overestimate the column density at (low 
S/N, low depth). 



8 



Fox, Savage, & Wakker 



FUSE/LiF; b=24 km b - '; depth = 0.50; S/N= 8; 500 runs 




Fig. 3. — Histogram of N a /Nt Tue for one model run, for three re cases. There is a significant dispersion in ovcrcstimation factors 

for a given case. As the rebinning increases, the most likely value of N a /Ntrue approaches unity 



